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We calculate the conductances of a Y-junction set-up of Luttinger liquid wires threaded by a 
magnetic flux, allowing for different interaction strength ^ g m the third wire. The scattering 
matrix and the matrix of conductances are parametrized by three variables. For these we derive 
coupled RG equations, up to second order in the interaction, within the scattering state formalism. 
The fixed point structure of these equations is analyzed in detail. For repulsive interaction [g, gs > 0) 
there is only one stable fixed point, corresponding to the complete separation of the wires. For 
attractive interaction {g < and/or ga < 0) four fixed points are found, whose stability depends 
on the interaction strength. For special values of the interaction parameters (a) g = 0, g^ < or 
(b) (?3 + g^ /2 — and g < we find whole lines of stable fixed points. We conjecture that lines of 
fixed points appear generically at interfaces in interaction coupling constant space separating two 
regions with different stable fixed points. We also find that in certain regions of the g-ps-plane the 
RG flow is towards a fixed point without chirality, implying that the effect of the magnetic fiux is 
completely screened. 



Introduction. Electron transport in strictly oncdimen- 
sional quantum wires is governed by the Coulomb in- 
teraction between electrons, found to destroy the Lan- 
dau quasiparticle concept valid in higher dimensions. An 
important part of that physics is captured by the ex- 
actly solvable Tomonaga-Luttinger liquid (TLL) model, 
in which backward scattering is neglected.'^ As is well 
known, transport through junctions of TLL wires is 
strongly affected by the interaction in the wires, some- 
times leading to a complete blocking of transmission. 
Many of the early works on the problem of two TLL- 
wires connecte d b y a junction employed the method 
of bosonizatiorPEi, or, for special values of the inter- 
action parameter, used a mapping on to exactly solv- 
able modelPI^. An alternative formulation using the 
fermionic representation has been pioneered bj^. There 
the renormalization group (RG) method has been used. 
The advantage of the latter formulation is that it is rather 
flexible, allowing for the inclusion of spin, of backward 
scattering, the study of w ires with two impuritie^, and 
of multi-lead junction^SHI in and out of equilibrium. The 
limitation of the original model calculation- to weak cou- 
pling has prevented a more widely spread use of this 
method. It may be shown, however, that the method 
may be extended to apply in the strong coupling regime 
as well, by summing up an infinite class of leading con- 
tributions in perturbation theorjffSl. Quite generally, the 
transport behavior at low energy/temperature is domi- 
nated by only a few fixed points of the RG flow. The 
most common and intuitively plausible fixed points are 
those with quantized conduction values, G ~ 0, 1, in 
units of the conductance quantum Go = e'^/h. There 
may, however, appear additional fixed points associated 
with noninteger conductance values. A particularly in- 



teresting case is that of a symmetric three-lead junction 
with broken chiral symmetry, as induced by a magnetic 
flux threading the junction. This latter proble m has been 
studied by the method of bosonization irfSM!^ and by the 
functional RG method in^^. These authors have identi- 
fied a number of new fixed points, in particular for at- 
tractive interaction. The discovery of these new fixed 
points has sparked interest in the problem of mapping 
out the complete fixed point structure of the theory, even 
though most of the interesting new physics appears to be 
in the physically less accessible domain of strongly at- 
tractive interactions. Only the fully symmetric situation 
o f iden tical wires {g = ga, see below) has been considered 
jjpUm^ As we show below, by restricting the consider- 
ation to the symmetric case, important new physics is 
missed. We mention in passing that a very general study 
of junctions of n quantum wires using bosonization and 
scale invariance can be found iiP. The effect of spin 
o n chir al Y-junctions has been considered very recently 

Experimentally the Y-junction set-up may be realized 
by a one-dimensional tunneling tip contacting a quantum 
wire. Although Y-junctions threaded by magnetic flux 
have not yet been studied in experiment (the flux gen- 
erated at a nanoscale junction by an externally applied 
magnetic field will be too small, but local magnetic mo- 
ments generating internal flux may be an option) , similar 
set-ups with either broken time reversal symmetry or bro- 
ken parity have been considered recently. The tunneling 
of spin-polarized electrons into the edge channel of a spin 
Hall effect system has been theoretically studied u^^. 
The asymmetry of the injected current calculated there 
is reminiscent of the "Hall current" (conductivity Gab) 
studied below. However, the effect considered irP3 in- 
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volves the symmetric component of the off-diagonal con- 
ductivity (in our notation Gab = Gba)- Similarly, by 
injecting current into the edge states of a semiconductor 
Hall bar in the presence of a suitable magnetic field, an 
asymmetry of the currents at the two half wires has been 
detected, depending on the interaction strength in the 
wir^i^. The latter has been interpreted as a signature of 
charge fractionalization. 

In this paper we report results of an RG treatment of 
electron transport in the linear response regime through 
a junction of three spinless TLL wires threaded by a mag- 
netic field. We employ a fermionic representation as de- 
scribed in detail irflSl. A comparison of this method in its 
simplest form with the functional renormalization group 
method has been made recently irP^. The method has 
formerly been used irP121 in the case of Y- junctions. More 
recently this approach has been used to derive the RG 
equations for the conductances of a Y-j unction connect- 
ing three TLL wires in the absence of magnetic flux, for 
weak interactiort^, and in a recent work for any strength 
of interactiorPi'. IrPSl it was found that even in weak cou- 
pling, but beyond lowest order, interesting new structures 
appear. The most important result there is that the fixed 
point found in lowest order for repulsive interaction, de- 
scribing the separation of the third wire (the tunneling 
tip) from the ideally conducting main wire, is actually a 
saddle point and is thus unstable. The message here is 
that even a weak higher order contribution may change 
the RG-flow in a dramatic way. In the present paper we 
describe a similar effect: an asymmetry of the three wires 
of a chiral Y-j unction to the effect that in the tip wire 
the interaction strength is different from that in the 
main wire, allows to access certain lines in interaction 
parameter space where a whole line of fixed points rather 
than a single fixed point is stable. This happens for at- 
tractive interaction only (g, g-^ < 0) The line of stable 
fixed points is connecting two fixed points at two special 
manifolds of interaction values, (a) g = 0; 53 < 0, or (b) 
whenever the condition (73 — ~g^ /2 is met. Along these 
fixed point lines the FP values of the conductances are 
continuously varying. To our knowledge this is the first 
time that lines of stable fixed points of the conductance 
have been found in models of TLL- wires. ^^Sl 

The model. We consider a model of interacting spin- 
less fermions describing three quantum wires connected 
at a single junction by tunneling amplitudes in the pres- 
ence of a magnetic fiux = ^/^q piercing the junction 
(with $0 the fiux quantum). In the continuum limit, lin- 
earizing the spectrum at the Fermi energy and including 
forward scattering interaction of strength gj in wire j , 
we may write the TLL Hamiltonian in the representation 
of incoming and outgoing waves as 

ifj"* = 2^TVFgJi>l,„iJJ^^ni^]^,^ti'J^out Q{x; -L, -I) . 
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FIG. 1; (Color online) (a) The geometry of a Y-junction 
is shown together with currents of incoming and outgoing 
fermions, the type of fermions in the scattered states repre- 
sentation is indicated by its color. The short-range interac- 
tion between the fermions is shown by wavy lines, (b) The 
equivalent representation in the chiral fermion basis after the 
unfolding procedure; the initially local interaction becomes 
non-local. 

Here Q{x; —L, —I) = 1, if —L < x < —I and zero other- 
wise, where L is the length of the interacting region in 
each wire and I is the size of the junction, inside which 
interaction is assumed to be absent. The regions x < —L 
are thus regarded as reservoirs or leads labeled j = 1, 2, 3. 
We assume that the boundaries at x = —L do not pro- 
duce additional potential scattering. The so-called back- 
ward scattering in the wires is unimportant in the con- 
sidered spin-less situation, as it can be included into the 
gj amplitudes.^ We denote gi = 92 — 9 from now on, and 
put vf — ^- The various incoming and outgoing channels 
are illustrated in Fig. [ij 

In terms of the triplet of incoming fermions ^' = 
("01,™! V'2,ira) '03,in) the Outgoing fermion operators may 
be expressed as '^{x) — S ■ ^(— x) . Here the 3x3 S- 
matrix characterizes the scattering at the junction and 
(up to irrelevant phase factors) has the structure 

/ ri, ti2, ii3\ 

S = t21, ri, t23 , 

\t3i h2 r2 j 
ri = i (cos6lH-e''^)sin^, r2 = cos6', (2) 
^13 =^32=* COS I sin 6^ ^23 = ^31 = * sin 6 sin | , 
ti2 = cos6'cos^ I - e^'^'sin^ |, ^21== ii2lc^7r-5 • 

We express the interaction term of the Hamiltonian in 
terms of density operators Pj^in — tpj ini^jM — Pj^ = 
Pj, and pj^out = '0],o«tV'j,o«t = = Pj, where pj = 

5+ • pj ■ S, as 

i/j"* = 2^{.g[pi^i + P2P2] + 93P3P3}Q{^; ~L, -I) . (3) 
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The density matrices are given by {Pj)ai3 = SapSaj 
and ipj)ai3 = S'^jSjp- A convenient representation 
of 3 X S-matrices is in terms of Gell-Mann matrices Aj, 
j = 0, 1, ... ,8, the generators of 5;7(3) (secP). Notice 
that the interaction operator only involves A3 and Ag (be- 
sides the unit operator (Ao)q,;3 ~ ^2/3(5^^). We note 
Tr{\j) = 0, Tr{\j\k) = 25jk. J = 1, . . . , 8, and [A3, Ag] = 
0. Introducing a compact notation we define a three- 
component vector A = (A3, Ag, Aq), in terms of which the 
densitries may be expressed as pj — •\/l/2^^ Rj^X^, 
where the 3x3 matrix R is defined as 
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and has the properties R "'^ = R"^, defR. = I. The 
outgoing amplitudes are expressed in the analogous form 
with \j replaced by \j = 5+ • Aj • 5 . With the aid 
of the \j the S-matrix may be parametrized by eight 
angular variables (se^^H) Yov the case under consid- 
eration only three of these, are relevant: S = 
^i\2mg^M^-4>)f,^\5e^^\2{n-0/2 . ^his S'-matrix is given 

explicitly by Eq. As will be shown below, the S- 

matrix and therefore the angular variables 9, ip, ^, will be 
renormalized by the interaction. 

Parametrization of conductance matrix. In the linear 
response regime, we may define a matrix of conductances 
Cjk relating the current Ij in lead j (flowing towards 
the junction) to the electrical potential Vk in lead k : 
Ij = J^k^jk^k- The conductance matrix C has only 
three independent elements, which may be represented 
in the form of a 2 x 2 antisymmetric matrix G with el- 
ements Gil = Gaa = |(1 - a), G22 = Gbb = |(1 - b), 
G12 = —G21 — Gab — -?=■ These conductances relate 

v3 

the currents la — I2), lb — ^{Ii + h — 2/3) to the 

bias voltages K = {Vi - F2), Vb = ^iVi + V2- 2V3). It 
is useful to note that the nonzero elements of the matrix 
— R"^ • C • R are essentially the reduced conduc- 
tances: Cfi = 2Gaa , G^2 = (3/2)G66 and C{^2 = ~Gii = 
\/3Gab , where the numerical factors arise due to the 
physically motivated asymmetric definitions of the cur- 
rents and voltages. We now show that the parametriza- 
tion of conductance in terms of a, 5, c follows quite natu- 
rally, by observing that the initial conductances are given 
by Gjk = Sjk - Tr{p^^pk) = 5jk - l-S'JfeP, where the la- 
bel r indicates that the quantity is fully renormalized by 
interactions.'^il By expressing Yjk = Tr{pjPk) in terms 

of Y^^ = lTr{y^X^) as Y = (R • • R^), we find by 
comparison with conductance matrix that Y^ has 
nonzero elements given by the conductance parameters 
introduced above V/l = a, = b, Y^ = -Y^ = c, 
¥33 = 1. Therefore, we may use Y^ to represent the con- 
ductances in the renormalization group analysis below. 
By substituting the S'-matrix in the form ([2]) into the def- 
inition of Y^, we find the relation of the a, 6, c to the an- 
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FIG. 2: (Color online) Allowed values of conductances, a, b, c, 
lie inside the depicted body. The location of RG fixed points 
N, A, X*, M is also shown. The line of fixed points connect- 
ing FPs A and x'^ is indicated in yellow. 

gle variables: a = —^{cos^ 6 + 1) cos^ ^ -t- cos cos ■0 sin^ ^, 

b = i(3cos2 6* - 1), c = sin^ 61 cos e We find there- 
fore that a, b, c are confined within the region a S [—1, 1], 
b e [-1/2, 1], c e [-V3/2, 73/2]. The physically allowed 
points in a-6-c-space are inside a complex-shaped body 
labelled B as shown in Fig. |2] 

Renormalization group equations. The renormaliza- 
tion of the conductances by the interaction is deter- 
mined by first calculating the correction terms in each 
order of perturbation theory. We are in particular inter- 
ested in the scale-dependent contributions proportional 
to A = ln{L/l) , where L and / are two above lengths, 
characterizing the interaction region in the wires. Here 
we only consider terms up to second order in 3,53. 

The conductances are given by Cjk = 

I-L'^yi(Pji-^) - PjiL))pk{y))uj^a- Here the com- 
bination pj{—L) — Pj{L) means that the current is 
measured at x — — L, and integration over y for the 
incoming densities corresponds to the notion of voltage, 
14, applied to the leads; in the considered limit ujL — )■ 
the value of L in this definition does not play a role.f^^ 
The quantum averaging (. . .) is accompanied by taking 
trace over the matrices' product. Up to second order we 
need to consider the contributions depicted in Fig. 4 of 
Ref.'22]. Differentiating these results with respect to A 
(and then putting A = 0) we find the RG equations up 
to second order in the interaction 

Her^the = (R^ • W • R)^^ = ^[A^, A,], and W^k = 
[pj, Pk], are 3x3 matrices for each pair of pu (or jk) and 
the trace operation Tr is defined with respect to that 
matrix space, whereas gj}^ = (R-^ • g • R)^^,, with gjk = 
Sjkgj and are scalars with respect to this space (sums 
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over repeated indices are implied). The nine matrices 
may be evaluated best with the aid of computer 
algebra, inserting the ^-matrix in terms of the quantities 
a, b, c, as given above. As a result one finds the following 
set of RG-equations 

^ - /3,= f(^2 + ii32 + fi?3), (5) 

If = /3, = -4[l + 2a+i(l + 26) + fC3], 
where we defined 

Al = f + fe-2a^ A2^a{l~b) + c^ , 
B2 = (l-fe)(l + 26), 

As ^ a{3 + b - Aa^) + + l{ab ~ a + {4a + 1)) 

+ ^{l + 2b){a-ab + c^), 
^3 = (1 - b) {2a^ + 2b+l) + 2ac^ 

+ l{2b+l) (6 + 2c2- 1) + Jj(i + 3fe-4&3), 
C3 = + 2a + 26 + f - f (a(l - 5) + 3 (6 + c^)) 

and K = 35/(5 + 253)- 

Fixed points of RG flow. The fixed points of the above 
RG equations may be found analytically. For general 
(but small) values of g, 53 we find five fixed points in 
total, labelled N, A, M, x+ , X~ • fixed point N we have 
a = b = 1 and c = 0. It corresponds to the complete 
separation of all wires and Gaa — Gtb — Gab — 0. The 
approach to the fixed point is given by the power laws 

N N 

Gbb oc L~'^^ and 4Gaa — Gbb oc L~^^ with exponents 
7f = 5+.93+.g^/2 and 7^ = 2g+g'^. The third direction 
of RG flow, along the c axis in Fig. [2j is unaccessible due 
to unitarity restrictions, se^^Sl^ This fixed point is the 
only stable one for repulsive interaction, either 5 > or 
53 > 0, or both. 

Fixed point A is characterized bya = — I, b = 1, c — 
corresponding to Gaa = 1, Gbb — Gab — and thus 
the situation of a perfect wire 1-2, detached from wire 
3. Stability analysis shows that it is stable for attractive 
interaction, g < 0, provided 53 + g^ /2 > 0. The critical 
exponents at A are defined as 4:{Gaa — 1) + Gbt oc L^'i'i 
and Gbb,Gcc oc L~^2 and are given by ■^f = 2g — g^, 
72^ = -53"5V2. ^ 

Fixed point M is characterized by noninteger values 
of the conductance, which are found to depend on the 
strength of interaction, a~ ~K/3+gh, 6~ (k^/3— l)/2— 
gnh, c — 0, with h = k{k — l)/6. It is stable for 5 > 
and 53 < in a region depicted in Fig. [3] The critical 
exponents at M are obtained as 'Ji^ ~ g{K — 9k^^)/12, 
-f^'^ ~ -g{K + 3k-1)/2 and 7^^ ~ 5(3 - k)/6, in the 
scahng laws ^Gaa + Ghb - (l + f )^ oc L"-^" , 4^Gaa + 

Gbb + « L-^" , and Gcc oc L'-'S' . The region 

of stability of M corresponds to < 5 < —53, which 
translates to —3 k < 0, in which range all 7*^ > 0. 



Whereas fixed points N, A, M exist also in the absence 
of magnetic flux, two new (stable) fixed points, x^i ap- 
pear as a result of the broken time reversal symmetry. 
These ar e the fixed points of maximum chirality first dis- 
cussed ir l^^l^^ l. The corresponding fixed point values are 
a = & = — i and c = ±^\/3 , or in terms of conductances: 
Gaa = f , Gbb = i, and Gab = -Gba = ±5- The fixed 
points X* are stable for attractive interaction, g, 53 < 0, 
in the domain I53I > 5^ (see Figj3|. In this case the 
conductances are not quantized in integer units of the 
conductance but assume fractional values. The two fixed 
points x"^ are located symmetrically to the plane c — 
in a-6-c-space. Depending on whether the initial condi- 
tion for c is Cin > or Ci„ < the RG-flow will tend 
to x+ or x~ ■ If — the flow will stay in the c = 
plane; the M point is stable in place of in this cas^^D, 
The critical exponents at x^ are found as 7J- = 5 and 
1^ =93+9^/2. 

In the case 93 — g the values of the critical exponents 
at the fixed points N, A, x* obtained here agree with 
the results quoted in the literatur d^^ l ^^ l The comparison 
of our results with those obtained by the bosonization 
method is nontrivial, since the latter are valid for infi- 
nite TLL wires, whereas our method provides results for 
the experimentally relevant case of finite TLL wires at- 
tached to ideal leads. Following earlier works, Oshikawa 
et aL^i^l have defined a prescription for calculating the 
conductances Gij of wires attached to ideal leads from 
the conductances Cij of the infinite wire system, see Sec. 
12 therein. In the following we shall use their recipe for 
converting Gij into C'ij. 

Fixed point M requires special c onsid eration. Its ex- 
istence had been first conjectured irfSE^l^ but its proper- 
ties were not accessible. Later fixed point M was clearly 
identified in numerical studies using the functional renor- 
malization group methocP^J. The above values of the crit- 
ical exponents [see alscPl^, 7^^ — —2g ~ 2{K — I) and 
73*^ = 5/3 ~ (1 - K)/3, with K = V(l-5)/(l + 5) the 
Luttinger parameter, are in agreement with the numer- 
ical estimates obtained irfi^. Recent numerical result^^D 
on the values of conductances Cij at x^ and M obtained 
via the Density Matrix Renormalization Group (DMRG) 
method confirm the known analytical results for the x^ 
fixed points, in agreement with our results. In Ref.'^il (see 
also above) we found the value of the conductances at the 
M point, Cij — |(3(5jj — 1), independent of the strength 
of interaction, g — gs- Expressing the results in terms 
of conductances Cij, we find Cij = Cij^^^. The com- 
ponent G12 agrees very well with the numerical estimate 
given in^^, lending support to our identification of the 
fixed point M with the "mysterious" fixed point irf^^l. 

It is worth noting that in the region of stability of the 
M point the initial Hall conductance scales to zero, mean- 
ing that the magnetic flux will be screened (this is also 
true when the or ^ points are stable, but for a trivial 
reason, since at least one of the wires gets separated in 
that case). In Fig. [2] we show the location of the fixed 
points in a-6-c-space for a typical choice of interaction 



FIG. 3: (Color online) The RG phase portrait, showing re- 
gions of stability of fixed points in the space of interactions g 
and gs. The M point only exists in the shaded regions. Lines 
of fixed points appear for g, 33 values located along the red 
(thick) lines. 



constants. All fixed points are located at the surface of 
the body B. 

Lines of fixed points. As discussed above, the line 
53 + 5^/2 — in the space of interaction constants sep- 
arates two domains in which fixed point x+ (or x ) 
or A are stable. For g, g^ on that line, in addition 
to the two fixed points A, we find a whole line of 
stable fixed points connecting these two. This means 
that the RG flow may stop at any point on this line, 
depending on the initial conditions. The fixed point 
line is parametrized by a = —1 -f ^r, h = 1 — |t, 

c = -^T, with T € [0, 1]. Substituting this parametriza- 
tion into the above RG-equations, we find ^{a, 6, c} — 
— 5(53 + g^/2)T(l — t){1, — 3, a/S}, demonstrating that 
on the line A-x^ £^nd for (53 + .g^/2) ^ the fiow is 
directed along the line, towards A or x^, depending on 
whether —173 < 5^/2 or the opposite. When g^ = —5^/2 
the RG fiow is directed perpendicular to this line with 
the exponents g{\ ± -^/l — t), i.e. the fixed points on the 
line are stable. This is demonstrated in Fig. [4] where 
the RG-trajectories on the surface of body B are shown 
for different g, g^ and different initial conditions. We ob- 
serve that the blue and green trajectories for interaction 
strengths in the regions of interaction parameter space la- 
beled x+ and A terminate at the respective fixed points 
(note that the line at 9 = n corresponds to a single point), 
whereas the red lines calculated for interaction values at 
the boundary between those two regions are found to 
terminate at any point along the lines connecting A and 

x+. 

There is a second line of fixed points, separating the 
domain of stability of the fixed points M [g > 0, 33 < 0) 
and x^'^ {9 < 0,53 < 0). It is realized for any 33 < 
at 5 = 0. The line of fixed points connects fixed points 
M and ^"^^ is parametrized by a = — ^t^, h = — ^, 



FIG. 4: (Color online) RG trajectories on the surface 7/) = 
red (solid): g — —0.333, (73 = —0.062, blue (dashed): g = 
-0.25, gs = -0.25, green (dotted): g = -0.4, 33 = -0.01. 



c = with T € [—1,1]. Finally, a trivial line of fixed 
points exists for any 33 > at g = 0. It connects the 
fixed points A and A'^ and corresponds to the interacting 
third wire, being detached from the non-interacting wire 
1-2 of arbitrary conductance, a G [—1, 1], 6 = 1, c = 0. 
Evidently, a should not be renormalized in this case. 

The present weak coupling theory may be extended 
into the strong coupling regime, using the approach 
confirming the existence of lines of fixed points. 

Conclusion. In this paper we developed a fermionic de- 
scription of a Y-junction (1-2 symmetric junction of three 
TLL-wires) threaded by magnetic flux, in the regime of 
weak interaction. In this case the conductance matrix 
has three independent components, the conductance as- 
sociated (i) with current flowing in wires 1 and 2, (ii) with 
current flowing from wire 3 into wires 1, 2, (iii) with cur- 
rent in wires 1, 2, driven by a voltage drop in wire 3 (Hall 
conductance). We calculated the scale dependent contri- 
butions to the conductances up to order g^ and 33 in the 
interaction constants. From there we derived three renor- 
malization group equations for the conductances. These 
equations have five fixed points labelled N, A, M, x"*", X~ ■ 
When both interactions are repulsive {g > and g^ > 0) 
only N is stable, which corresponds to the total separa- 
tion of the wires. If at least one interaction is attractive 
{g and/or g^ < 0), we find four regions in g^-g-space in 
which N, A,x^,M are stable, respectively. At the two 
border-lines separating regions of stability of A, x^ (the 
line gz+g'^/2 = 0) and of x'^,M (the line 5 = 0, 53 < 0), 
respectively, we find a novel and unexpected situation: 
at these lines in g^~g-s])ace not only the two stable fixed 
points assigned to the adjacent regions are stable, but 
a whole line of fixed points connecting those two fixed 
points becomes stable. We have checked that the leading 
higher order terms do not destroy the fixed line prop- 
erty (to be published). We conjecture that this is an 
incidence of a more general phenomenon: whenever two 
regions in interaction constant space with different sta- 
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ble fixed points touch, and the two fixed points neither 
vanish nor merge at the contact manifold, a whole line of 
stable fixed points connecting the two initial fixed points 
emerges. 
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